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. Abstract. Using the ~ 6 pb 1 of e + e annihilation data taken at t/)'(3686) with CLEO III 

O ' and CLEO-c detectors with estimated ~ 3.0 x 10 6 tp' events, we have searched for the /i c (l 1 Pi) 

| state of charmonium in the reaction ^'(3686) — > n°h c — > (77) (777c). The preliminary results are 

■ reported. 
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1. Introduction 

Charmonium spectroscopy has played a crucial role in the understanding of the quark-gluon 
structure of hadrons and the underlying theory of Quantum Chronodynamics (QCD). This is 
t— I . primarily due to the fact that the charmonium system is expected to be far less sensitive to the 

problems associated with relativistic effects and the large value of the strong coupling constant, 
a s , than the light quark (u,d,s) systems. Formation cross-sections for charmonium states, their 
masses and widths are also favorable for precision measurements. The existing experimental 
data have defined the spin-independent one-gluon exchange part of the qq interaction quite well, 
however, the spin dependence of the qq potential is not very well understood. In particular, 
the si • s~2 spin-spin, or hyperfine interaction is not well understood, because there is little 
experimental data to provide the required constraints for theory. The primary experimental 
data required for understanding the qq hyperfine interaction is hyperfine, or spin-singlet/spin- 
triplet splitting: AM hf (nL) = (M(n 3 Lj)) - M(n l L J=L ). 

For nearly 20 years, the only hyperfine splitting known was that for the IS" states of 
charmonium, AM hf (lS) = M(J/f) - M(rj c ) = 116 ± 2 MeV. Very recently, Belle, CLEO 
and BaBar succeeded in identifying rj' c (2S), with the rather surprising result that AM] % f{2S) = 
M{ip') — M(ji' c ) = 48 ± 4 MeV. Potential model and quenched lattice calculations predicted a 
larger AM hf (2S) [1]. 

It is of great importance to find out how the hyperfine interaction manifests itself in P states, 

1. e., to find AM hf (lP) = M(< 3 Pj >) - M( 1 Pi). With scalar confinement, AM hf (lP) = is 
expected. It is necessary to determine if this is true. The c.o.g. of 3 P states, M(< 3 Pj >), 
is well measured, M(< 3 Pj >)=3525.3±0.1 MeV. What is needed is to identify h c and make a 
precision measurement of its mass. 

2. Prior Experimental Searches for h c 

The Crystal Ball experiment at SLAC made a search for h c in 1982 [2]. The search was 
unsuccessful and they reported 95% confidence limits of B(ip' — > ir h c , h c — > jr] c ) < 0.32% 
in the range Mh c = 3440 — 3543 MeV. The next search for h c was made by the Fermilab 
experiment E760 [3] in the reaction pp — > h c — > ir°J/t(j. It was claimed that a statistically 



significant enhancement was observed and that the data indicated M(h c ) = 3526.2 ± 0.15 ± 0.2 
MeV. However, such an enhancement has not been confirmed by the successor Fermilab E835 
experiment, with significantly higher statistics [4,5]. The E835 experiment also searched for h c 
in the reaction pp — > h c — > ^r] c . Preliminary evidence at the ~3<7 significance level has been 
recently reported with M{h c ) = 3525.8 ± 0.2 ± 0.2 MeV [5]. No positive evidence has been 
reported yet by Belle and BaBar Collaborations. 

It is fair to say that at present there is no convincing experimental evidence for h c observation. 

3. CLEO Searches and Results 

The above considerations have motivated us to search for h c in the ~ 6p6 _1 data taken at CLEO 
with estimated ~ 3.0 x 10 6 ip' events, in the reaction 

y/ -> TT°h c , h c ->■ 7%. 

We search for this channel: (a) without using r\ c decays (INCLUSIVE approach, see Section 
3.1), and (b) using six dominant r/ c decay modes (EXCLUSIVE approach, see Section 3.2). In 
both methods we search for h c in the mass recoiling against ir° from decay ip' — > ir°h c . This 
method benefits from the excellent resolution of the CLEO calorimeter. 

3.1. Inclusive Analyses 

Two independent analyses have been performed, and results from the two are consistent. I will 
describe one of them in detail, and will later mention the differences between the two analyses. 
We use the following selection criteria: N s hower >3, N tra ck >2. The selection of the showers and 
charged particles are done using the standard CLEO quality cuts. 

We reconstruct 7r°'s by requiring that the two photon invariant mass be in the range 
M 77 =135±15 MeV, and that the two photons have been succesfully fitted to tt°. We require 
that there be only one 7r° in the event with a recoil mass in the expected h c mass range of 
3526±30 MeV. 

The ip' — > 7r + 7r~ J/ip and ip' — > ir°ir°J/i() events are removed by cutting on the recoil mass of 
7r + 7r~ and 7r°7r°, respectively. 

We define hard 7's, the possible candidates from h c — > 7^ c decays, by > 400 MeV. We 
reject such 7's which make a 7r° or r\ with any other 7's. We then require that the energy of 
hard 7 should be in the range .E 7 =503±40 MeV. 

The background in data has been fitted in three ways: (a) ARGUS shape, y = x x 
y/1 — (x/a) 2 x exp(b x [1 — (x/a) 2 ]), (b) second-order polynomial shape, (c) background shape 
from Monte Carlo. The significance levels are obtained as a = y/— 2 \n(Lo / L max ) , where L max 
and Lq are the likelihoods of the fits with and without the h c resonance. 

The analysis on the Monte Carlo samples has been performed. The event selection criteria 
applied to the Monte Carlo samples were identical to those applied to the data. 10,000 signal 
Monte Carlo events for the channel tp' —* 71-0 h c — ► (77) (iVc) were simulated. The recoil mass 
distribution against 7r° in signal Monte Carlo, for input T(h c )=0 MeV is well fitted with a double 
Gaussian with parameters oi=1.3 MeV, CJ2=3.7 MeV, and the fraction of second Gaussian was 
0.43. These parameters, which represent the 7r° recoil mass resolution at h c , are used to fit 
the signal in the data. The selection efficiency was about 16%. We also analyzed a sample of 
~ 12 x 10 6 generic tp' Monte Carlo events (events containing all measured ij)' decays except those 
via h c ) in four separate samples, each with approximately the same size (~ 3 x 10 6 ) as the data. 
The signal Monte Carlo events were added in to the generic Monte Carlo. The study of these 
Monte Carlo events yielded good agreement between input and output values for both, M{h c ) 
and B(tp' — > n h c ) x B(h c — > 7%), and showed that the analysis is sensitive to h c production. 

Figure 1 shows recoil mass distribution against tt° in data. The results of the fit are: 
M(h c ) = 3524.4 ± 0.7 MeV, N(h c ) = 156 ± 48, significance (h c ) = 3.3 a. 



An independent alternative analysis has been done. The main difference is that in this 
analysis instead of constraining the energy of the hard photon, the constraint is put in terms 
of recoil against 7r° 7 (jic mass). The results are consistent with those shown above. Thus our 
preliminary CLEO results from two inclusive analyses are: 

• M(/i c )=3524.8±0.7(stat)± ~l(syst) MeV, 

• Btf TT°h c ) x B(h c -> 7%) =(2-6) xl(T 4 , 

• The significance of h c detection > 3 a. 

Estimates of systematic errors in M(h c ) have been made by studying the following: tt° energy 
scale, background shapes, Monte Carlo input/output differences, non-resonant background, 
assumed h c width, binning effects, cut variations, and finally, the difference in M{h c ) in the 
two inclusive analyses. 

\\i / ^%°h c ^(yY)('m c ) 



Figure 1. Distribution of the recoiling mass 
against tt° in data (inclusive analysis). The 
curves are the results of the fit. The shape 
3.5 3.5i 3.52 3.53 3.54 3.55 of the signal is assumed as Double Gaussian, 

7i° recoil mass in GeV and the shape of the background is assumed 

as ARGUS shape (see text). 



3.2. Exclusive Analysis 

Six rj c decay modes which have reasonably high PDG04 branching ratios have been studied: 
KgK^iT^, K + K~ir°, K + K~ir + 7r~ , 27r + 27r~, 7r + 7r~77 (77 — > 77), and 7r + 7r — 77 (rj — > 7r + 7r~7r ). 

Standard CLEO selections are used for showers, tracks, and particle identification. The total 
energy-momentum conservation of the event has been required, and the invariant mass of the 
77 c decay candidates are required to be close to the nominal r? c mass (within 50 MeV). Figure 
2(upper plot) shows the ir° recoil mass distribution for the sum of the six exclusive channels. 
The fit results are: 

• M(/i c )=3524.4±0.9(stat) MeV, 

• N(/i c )=15.0±4.2, 

• The significance of h c detection ~ 5 a. 

Note that the significance is calculated using likelihood differences. The background estimation 
by using rj c sidebands(closed circles in Figure 2, lower plot), or by using generic Monte Carlo 
events(open squares in Figure 2, lower plot), yield consistent results. No estimate of the 
systematic uncertainty in M(h c ) has been made so far. 
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Figure 2. Distribution of the recoiling mass against 7r° (exclusive analysis). 



4. Summary 

We have analyzed ~ 3.0 x 10 6 ip' from CLEO III and CLEO-c to search for h c ( l Pi) production 
in the reaction tj/ — > TT°h c , h c — > ^r\ c by two methods 1. INCLUSIVE - which does not use r\ c 
decay modes, 2. EXCLUSIVE - which uses six hadronic decay modes of 7/ c . 
In the recoil mass spectrum of 7r°, we see an enhancement in both analyses. 

• In the inclusive analysis we obtain 

M(/i c )=3524.8±0.7(stat)± ~l(syst) MeV, 
B(i// -> 7r h c ) x B(h c 7 r/ c ) =(2-6) xlO- 4 , 
significance of h c detection >3 a. 

Thus, AM hf = (M(xj)) - M( x Pi) = 0.5±0.7(stat)± ~l(syst) MeV. 

• In the exclusive analysis we obtain 

M(/i c )=3524.4±0.9(stat) MeV, 
significance of h c detection ~ 5 a. 

• The inclusive and exclusive results for M(h c ) are in excellent agreement. 
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